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Operación de grupo aditiva + 

Se define la operación matricial + como A + B=A/BllB/A, dadas A y B matrices cuadradas de 
mismo rango con valores booleanos y, por tanto, A + B equivale a hacer una operación XOR 
elemento a elemento. 

Ejemplo: 



Observamos que la operación aditiva puede servir para comprobar si dos matrices son iguales, 
sin ir más lejos observamos que es fácil determinar su elemento neutro: 


Para cualquier A dentro del grupo, e + A = A + e = A, donde e = (o 1 1; ,} un o en todas las 
posiciones de la matriz. 


Ejemplo para matrices de rango 3 : 



Dada esta definición, el elemento opuesto de cada matriz es la propia matriz. 


A + (-A) = (-A) + A = e 
A = -A 


Así queda definido el grupo aditivo + para matrices booleanas cuadradas de mismo rango. 



Operación monoide multiplicativa x 


Definimos la operación de grupo multiplicativa entre dos matrices booleanas cuadradas de 
mismo rango de la siguiente manera: A x B = A & ((A-B)>2) 


Esta operación se compone de tres operaciones: 

i) A-B es el producto matricial convencional a través de productos escalares de filas por 
columnas. Transforma dos matrices booleanas en una matriz de enteros. 

Ejemplo: 




m-ci* ^ 
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2) XA: A > 2 es una función que transforma una matriz de enteros en una matriz de booleanos, 
sabiendo que el elemento en la posición i, j valdrá i si éste es mayor que 2. 

Ejemplo: 



> 2 



3) A & B es la función que calcula AND elemento a elemento entre dos matrices: el producto 
booleano no matricial. 

Ejemplo: 


1 

1 



Podemos encontrar el elemento neutro de esta operación siempre y cuando las matrices 
cuadradas de mismo rango y valores booleanos tenga un formato específico; 



Definición de matriz lógica: 

1) La matriz se compondrá de matrices booleanas donde en la diagonal principal habrá matrices 
identidad de diversos rangos. 

2) Los menores i, j que se ubiquen en la matriz triangular superior (i<j) serán matrices 
booleanas transpuestas a los menores j, i ubicados en la matriz triangular inferior. Y serán 
distintos de la matriz o. 

Este conjunto de matrices poseen como invariante una fórmula booleana, como más adelante se 
detallará. 


Ejemplo: 
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El elemento neutro para este conjunto de matrices 
sus menores. 

Ejemplo: 

/ 1 1 
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es el que tiene 1 en todas las componentes de 
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Para este elemento neutro no hay elemento opuesto, por lo que no buscamos que nuestro 
monoide sea un grupo. Observamos que, efectivamente, dependiendo del formato de menores 
dentro de nuestra matriz se conformará un monoide u otro. 



Monoide x de matrices lógicas 

Para que el operador multiplicativo sea un monoide debe: 

1) Cumplir la propiedad asociativa 

2) Tener un elemento neutro 

El elemento neutro ya fue mostrado en la sección anterior, sin embargo, dentro del invariante de 
las matrices, cada elemento ejerce sobre sí mismo de elemento neutro: A x A = A 

Por lo que pasamos a la cuestión de la asociatividad. Cuando las matrices no poseen un mismo 
invariante la operación no es asociativa; por lo que debe de aplicarse según la sucesión 
ascendente: 1-M, 2-M = i-M X 1-M, 3-M = 2-M X 1-M,... 

En lo que se refiere a la operación n-M x m-M = (n+m)-M se comprueba que es asociativa por la 
sencilla razón de que cada operación que se haga con la matriz M dentro de la sucesión genera 
un nuevo término en la sucesión como operación sucesor, por lo que la asociatividad se ve 
explicada a través de la asociatividad de la suma. 

Operación distributiva entre el monoide x y la suma + de matrices lógicas 

Vamos a comprobar cómo se da la operación distributiva entre tres matrices lógicas cuando se 
operan: Ax(B+C) = AxB + AxC.De esta manera quedará certificado que estos dos 
operadores conforman con tales matrices un anillo. 

La propiedad distributiva, no se dará para cualquier matriz, sino exclusivamente en las matrices 
que comparten el mismo invariante lógico. Por eso mismo habrá que usar la definición: 

A X B = n-M X m-M = (n+m)-M 

Asimismo, recordamos que A x B = A & ((A-B)>2), que significa determina la triangulación de B 
que rechaza casos en A. Por lo que A x (B+C) significa: determina la triangulación de lo que 
diferencia B con C para rechazar los casos de A, que sería como establecer las diferencias entre el 
establecimiento de la triangulación de B sobre A y de C sobre A:AxB+AxC. 

B + C = b-M + c-M = b-M / c-M cuando b < c, son los unos que están en B y que se dedujo que 
eran incoherentes con el invariante al madurar hasta C. 

Proposición. A x (B+C) = AxB + AxC para matrices de mismo invariante lógico. 
Demostración. 

La ausencia de esta demostración no desmerece las posteriores. De hecho, a pesar de que no 
tengo demostrado que se conforma un anillo, en el resto del documento se hablará como si lo 
fuera - considerando que también ha quedado corroborado en mis simulaciones el que se cumple 
la propiedad. 



Sucesiones ascendentes antes de llegar al ideal 

Mediante una matriz inicial M, cualquier multiplicación sobre sí misma nos dará el sucesor: 

M, 2-M = MxM, 3-M = 2 -MxM, ... 

Definimos la relación entre matrices lógicas: AdB cuando sum(A)> sum(B), esto es, cuando la 
suma de unos de la primera sea mayor o superior a la segunda. 

Proposición. ki-M 3 k2-M sii k2 > ki. 

Demostración. 

Partimos del coeficiente ki = 1 con M la relación entre las cláusulas dos a dos. A medida que se 
aplique la operación A i+1 = A¡ & ((A¡-B)>2), observamos que el número de í's va decreciendo y 
queda, por tanto, contenido en coeficientes mayores menores cantidades de í's. 

c.q.d. 

Cuando ki-M = k2-M con k2 > ki, entonces ki-M es el ideal. 

Desde la primera matriz M hasta la última antes de dar con el ideal se debe dar en todo el 
proceso: Bi d B2 d B3 d ... 3 Bn = Bn+i 

Proposición. Si k-A x A = k-A entonces k-A es el ideal. 

Demostración 

El algoritmo A i+1 = A¡ & ((A¡-B)>2) genera una secuencia Ai 3 A2 3 A3 3 ... 3 An donde el 
siguiente tendrá menos unos. Cuando en algún i A i+1 = A¡ entonces observamos que se repetirá 
exactamente la misma fórmula y, por tanto, entraremos en la condición del ideal, donde 
A i+1 + A¡ = o, que es el neutro de la operación aditiva. 

c.q.d. 

Proposición. La máxima longitud de una sucesión ascendente para matrices de rango n será 
siempre inferior a n 2 . 

Demostración. 

Dada la secuencia Ai 3 A2 3 A3 3 ... 3 Am, la secuencia más larga deberá eliminar, al menos 
un uno. Como el rango de la matriz es n eso quiere decir que no podrá eliminarse ni siquiera n 2 
considerando que las matrices en la diagonal principal poseen varios ceros que siempre quedan 
intactos. 

c.q.d. 

Para determinar el ideal sólo hay que dar con una posición lo suficientemente alejada del 1, 
esto es: M, 2-M= M X M, 4-M= 2-M X 2-M, 8-M= 4-M X 4-M,... 

A esta sucesión la llamaremos sucesión ascendente mediante duplicados 


Proposición. El número de pasos necesarios para alcanzar el ideal mediante duplicados siempre 
será menor de 2-Llog 2 nJ. Donde n es el rango de la matriz lógica. 



Demostración 

De la demostración anterior, se deduce que si el algoritmo busca el ideal mediante 
la secuencia M, 2-M, 4-M, 8-M,... entonces alcanzarán su destino siempre antes que una 
proporción logarítmica de veces la longitud más larga de la secuencia, que es n 2 , por tanto se 
deduce 2-Llog 2 n_|. 

c.q.d. 

Considerando que el invariante de la matriz de rango n tiene un tamaño de entrada para una 
máquina de Turing proporcional al rango de la matriz, queda demostrado: 

Proposición. Es necesario un tiempo de cómputo de 0 (x 2 log 2 x) para determinar si una fórmula 
booleana se satisface con x el tamaño de la fórmula escrita en la entrada de una Máquina de 
Turing. 

Demostración. 

De la proposición anterior, partimos de la representación de una fórmula bien formada como 
producto de alternancias, creando las variables temporales que se consideren oportunas. Como 
el número de variables no influye, al final consideraremos el número de cláusulas, que será 
proporcional al tamaño de la entrada y, a su misma vez, proporcional al rango de la matriz que 
representa el producto de alternancias. Por lo que cada operación matricial de rango x 
multiplicado por el número de veces que debe ser realizada la operación justifica 0(x 2 log 2 x). 
c.q.d. 

Mecanismo para generar la matriz inicial. 









@staticmethod 

def cM(clauseC, clauseR): 

R = matrix([[i] * len(clauseC)] * len(clauseR)) 
for i, X in enumerate(clauseC): 
if X in clauseR: 
j = clauseR. index(X) 

R &= matrix([[int((i==col) = (j==row)) \ 
for col in range(len(clauseC))] \ 
for row in range(len(clauseR))]) 
if -X in clauseR: 
j = clauseR. index(-X) 

R&= matrix([[int((i==col) A (j = =row)) \ 
for col in range(len(clauseC))] \ 
for row in range(len(clauseR))]) 

return R 


@staticmethod 
def SM(cl): 

M = [] 
for I in el: 

F = [[] for X in range(len(I))] 
for J in el: 

Ideal.addFilas(F, Ideal.cM(J, I).tolist()) 
M.extend(F) 
return matrix(M) 

@staticmethod 
def addFilas(F, M): 

'Prolonga en F las filas con M' 
for i in range(len(F)): 

F[i].extend(M[i]) 


Para inicializar: 


class Ideal: 


duplo = lambda M: M & (M * M > 2) 

def clone(self): 


return Ideal (clausulas = self.cl, 

def_init_(self, clausulas = [], M = None, A = [], I = [], denso = True): 

M = self.M.copyO, 

self.M = M if M is not None else Ideal.SM(clausulas) 

A = self.A.copy(), 

self.cl = clausulas 

I = self.I.copyO, 


denso = self.denso) 

self.card = [0] 


for X in clausulas: 

def_repr_(self) : 

self.card.append(self.card[-i] + len(X)) 

if not self.denso: 


return repr(self.cl) + "\n" + repr(self.M) \ 

self.V = self. variables () 

+ "\n" + repr(self.S) 


else: 

if not A and not I: 

return repr(self.cl) + "\n" + repr(self.M) \ 

self. A = set([]) 

+ "\n" + self.formatoQ 

self.I = set([]) 


self._actualiza() 


self._evalua() 


else: 


self.A = A # activos 


self.I = I # inactivos 


self.S = self._actualizaSalida() 


self.denso = denso 



Consulta y asignación de variables: 


def_getitem_(self, v): 

def_setitem_(self, v, valor): 

el, lit = self.Vfv] 

el, lit = self.V[v] 

M = self.menor(cl, el) 

self.selecciona(cl, lit, 

if M [lit, lit] == 0: 

activo = (valor= = i) A (self. el [el] [lit] < 0), 

return 2 if self.el [el] [lit] > 0 else 1 

actualiza = True) 

if M.sum() == 1: 


return 1 if self. el [el] [lit] > 0 else 2 


return 3 





Tratamiento sobre los literales: 



def selecciona(self, clausura, literal, 

def _actualiza(self): 

activo = True, actualiza = False): 

K = i 

'Pone a i el literal de la clausura de el en M' 

while K < 2*len(self.cl): 

if activo: 

K*= 2 

for i in range(literal): 

self. menor (clausura, clausura) [i, i] = o 

self.M = Ideal.duplo(self.M) 

for i in range(literal+i, len(self.cl[clausura])): 

def _actualizaSalida(self): 

self. menor (clausura, clausura) [i, i] = o 

self._evalua() 

else: 

self.S = dict([(V, i) for V in self.A if V not in self.I] \ 

self.menor(clausura, clausura) [literal, literal] = o 

+ [(V, 2 if V not in self.A else o) for V in self.I] \ 

if actualiza: 

+ [(V, 3) for V in self.V.keys() \ 

self._actualiza() 

self._actualizaSalida() 

if V not in self.A and V not in self.I]) 

def _evalua(self): 

def menor(self, i, j): 

'Comprueba sobre la matriz si las variables están activadas' 

return self.M[self.card[i]: self.card[i+i], self.card[j]: self.card[j+i]] 

'PRE: debe haberse actualizado la matriz' 

'POST: actualiza self.A y self.I' 
for v in self.V.keys(): 

Z = self[v] 

if Z = = i: 

self.A.add(v) 

elif Z = = 2: 

self.I.add(v) 





























Observa un caso: 



def_cali_(self, n): 

def formato(self, *rango): 

'Genera una observación posible' 

if not rango: 

other = self.cloneO 

rango = range(i, i + max(self.V.keys())) 

R = set([]) 

formatea = True 

Fi = o 

else: 

for i in range(len(other.cl)): 

formatea = False 

Mi = other. menor (i, i) 

self._actualizaSalida() 

if Mi.sum()==i: 

R = " 

X, Y = Mi.nonzeroO 

for K in rango: 

R.add(other.cl[i][X[o]]) 

R += repr(self.S.get(K, 3)) 

R.add(other.cl[i] [Y [o] ]) 

if formatea and K % 5 == 0: 

continué 

R += '•' if not K % 10 = = 0 else ' ' 

X, Y = Mi.nonzeroO 

N = len(X) 

X = X[n o/o N] 

Y = Y[n o/o N] 

R.add(other.cl[i] [X]) 

R.add(other.cl[i] [Y]) 

n //= N 

for z in range(len(other.cl[i])): 

ifz == X: 

continué 

other.M[Fi + z, :] &= o 
for z in range(len(other.cl[i])): 

if z == Y: 

continué 

other.M[:, Fi + z] &= o 

other._actualiza() 

Fi += len(other.cl[i]) 

return other 

return R 



































La máquina que edita el invariante: 


El constructor 


class Invariante: 

def_repr_(self): 

'Edita el invariante' 

return "Invariante: "+repr(self.F)+"\n"+repr(self.F) 

'0 siempre vale 3' 


' 1 siempre vale 1' 

def_getitem_(self, pos): 

' El opuesto de la posición es el valor opuesto' 

return opuesto(self.A.get(-pos, 3)) if pos < 0 else self.A.get(pos, 3) 

def_init_(self, nVariablesMax): 


self.t = nVariablesMax + 1 

def_setitem_(self, pos, val): 

self.AND = {} 

if abs(pos) > 0: 

self.XOR = {} 

val = val if pos > 0 else opuesto (val) 

self.F = [] 

self.A[abs(pos)] = self[pos] & val 

self .A = {1: 1} 

if self [pos] == 0: 


raise SinSolucion 

def clona(self): 

self.reduceO 

other = Invariante(self.t -1) 


other.AND = self.AND.copyO 

class SinSolucion (Exception): 

other.XOR = self.XOR. copy() 

pass 

other.F = [X[:] for X in self.F] 


other.A = self.A.copy() 

opuesto = lambda N: N//2 + 2*(N%2) 

return other 



Cómo se reduce la fórmula de manera superficial 


def _reduce(self): 

'A partir de las asignaciones sobre las variables se reduce la formula' 
'Devuelve las nuevas asignaciones, que no han repercutido' 

N = {} 
aElimF = [] 

for c, choice in enumérate(self.F): 
aElim = [] 

for 1, literal in enumérate (choice): 
if self [literal] == i: 

for 12 in range(len(choice)): 
if I2 == 1: 
continué 

N[abs(choice[l2])] = 2* int(choice[l2] > o) + int(choice[l2] < o) 
aElim.append(l2) 
break 

elif self [literal] ==2: 
aElim. append(l) 
if aElim: 
while aElim: 

choice.pop(aElim.popO) 
if len(choice) > 1: 

F[c] = choice 
if not choice: 

raise SinSolucion 
elif len(choice) == 1: 
if self [choice [o]] ==2: 

raise SinSolucion 
elif self [choice [o]] ==3: 

N[abs(choice[o])]= 2* int(choice[o] < o) + int(choice[o] > o) 
aElimF.append(c) 
while aElimF: 

self.F.pop(aElimF.popO) 
return N 


def reduce (self): 
N = self.A 


while N: 

N = self._reduce() 
for X in N.keys(): 
self.A [X] = N[X] 




Representación de todas las fórmulas 


def _temp (self, N = i): 

def opAND(self, Xi): 

self.t += N 

'Xi es una lista' 

return tuple(self.t - X for X in range(i, N + i)) if N > i else self.t -1 

Xi.sort() 


R = Xi.popO 

def sumador(self, X, Y): 

for X in Xi: 

if X > Y: 

R = self.sumador(X, R)[i] 

return self. sumador (Y, X) 

return R 

if X = = Y: 


return -i, X 

def opOR(self, Xi): 

if self[X] == i: 

'Xi es una lista' 

return -Y, Y 

return -self.opAND([-X for X in Xi]) 

if self[X] ==2: 


return Y, -1 

def opXOR(self, *Xi): 

if self [Y] == 1: 

'Xi es una lista' 

return -X, X 

Xi.sort() 

if self [Y] ==2: 

R = Xi.popO 

return X, -1 

for X in Xi: 

XaY = self.AND.get((X, Y), self._temp()) 

R = self.sumador(X, R)[o] 

XxY = self.XOR.get((X, Y), self._temp()) 

return R 

if not (X, Y) in self.AND.keys() or not (X, Y) in self.XOR.keys(): 


Zi, Z2 = self._temp(2) 

def SAT(self, * clausulas): 

self.F.extend([[-X, Zi, XaY], [XaY, Z2, -Y], [Zi, Z2, -XxY]]) 

'Cada argumento es una lista a satisfacer' 

self.AND[(X, Y)] = XaY 

for el in clausulas: 

self.XOR[(X, Y)] - XxY 

self[self.opOR(cl)] = 1 

return XxY, XaY 



Cálculos de vínculo aritmético 


def menor(self, N, M): 

def atribCostes(self, costes, n = None): 

’sum(N[i]*2**i) < sum(M[i]*2**i)' 

'La asignación de costes es a variables {V: C}' 

L= □ 

cOutf = lambda L: [ANDf(L[j], L[i]) for i in range(i, len(L)) for j in range(i)] 

Eq = 1 

ANDf = lambda Li, L2: L2 if Li == (1,) else Li if L2 == (1,) else tuple(set(Li + L2)) 

for i in range(len(N) -1, -1, -1): 

XORf = lambda Li, L2: tuple(set(Li + L2) - ({(1,)} if (1,) in Li and (1,) in L2 else set({}))) 

A, B = self.sumador(-N[i], M[i]) 
L.append(self.opAND([B, Eq])) 

cifra = lambda X, n: bin(X)[-(n+i)] == '1' if n< len(bin(X)) - 2 else False 

Eq = self.opAND([A, Eq]) 

if n is None: 

return self.opOR(L) 

n = 1 + int(log(supremo, 2)) 

N = [] 

for i in range(n): 

N = XORf(cOutf(N), [(V,) for V in costes.keys() if cifra(costes[V], i)]) 
return list(costes.keysO), N 



